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Goals & Key Features of TheSIS 
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GOALS 

1) Prepare for multivariate modem control loop design using detailed models of 
device and system dynamics, i.e. “thermal signatures”. 

2) Capture the maximum thermal performance information in the minimum time. 

KEY FEATURES 

1) Measure the thermal signatures of optoelectronic components to 
direct sequence spread-spectmm temperature cycles, 

2) Model the thermal signatures, and 

3) Select the most likely model per the Akaike Information Criterion (AIC). 

4) Analyze parameter variability and develop multivariate Figures of Merit. 
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Using the AlC-tested model and parameter vectors from TheSIS, one can analyze changes in model 
parameters, detect subtle reversible shifts in performance, investigate the cause of instabilities or 
irreversible changes in component or subsystem performance, e.g. aging, and select high-performing 
components on a multivariate basis, i.e. with multivariate Figures of Merit (FOM). We show examples 
of the TheSIS methodology for passive and active components and systems, e.g. fiber Bragg gratings 
(FBGs) and DFB lasers with coupled temperature control loops. 
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Applications 
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• EO Component, Subsystem, and System Evaluation (this talk) 

• Electronic Component Evaluation 

• e.g. Oscillators, Crystals, Frequency References 

• Reliability/Failure Analysis, 

• e.g. modify Thermal Vacuum (TVAC) cycling profiles to include System 
Identification (SI) via direct sequence spread spectmm “dither”. 

• Efficient Thermal cycling and Control Loop SI Problems 

• e.g. Tuning PID or other control loops -> “plant” models 

• e.g. Polymerase Chain Reaction (PCR) control 



Three Perfect Barker Codes 

Use for : Short Direct Sequence Spread Spectrum 



Goddard Space 
Flight Center 


Barker 3 : 
Barker 7 : 
Barker 11: 


+ + - 

+ + + ..I+ . 

+ + + + . _ + _ 


Each + or - is a “Chip” 



Note: the spectrum of the other perfect code, Barker 2 +-, is not shown 



Modified Barker 7 Profile 

(7 sub-intervals = "Chips") 
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Original -(4 6 7) Modified =(3.925 5.875 7) 



ChipXumber n 



Autocorrelation of Barker 7 
and Modified Barker 7 
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Original = (4 6 7) 



Modified =(3.925 5.875 7) 



Application: Active Optoelectronic System 
DFB Environmental Perturbation, T mount 
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GHz Frequency 5400 samples 



Note Frequency is from Laser Heterodyne method with the reference laser or LO at a 
higher frequency than the device under test (DUT). Therefore, as displayed here, a 
frequency increase corresponds to a redshift of the DFB laser. 


T chip= 1 minute 
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A Comparison of Autoregressive Moving Average (ARMA) Models 
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of DFB Laser Dynamics using the Akaike Information Criterion (AIC) 






e 


( AIC AlCmin ) /2 


is the likelihood of a 
given model relative to 
the most likely model. 
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Pole and Zeros of ( DFB Laser + Laser TEC) 
Response to Mount TEC set point changes 
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(AR1MA8; one of the seven Zeros is not shown; (at -3.48) ) 

Sampling Period, T=4.5 seconds 



-1.1 1.1 
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1. False Hysteresis 

2. Recursive peak finding algorithm 

(X c , Insertion Loss (dB), Reciprocal Bandwidth Squared) 

3. Modified Barker 7 Cycling 

4. FBG AR3MA2 Poles and “thermality zero” for MA2 models 

5. Stability of Poles and “thermality zero”, Z T 

6. Temperature coefficient of center frequency, a. 

7. Contrast conventional vs TheSIS temperature cycling 

8. Multivariate selection metric 



Center Wavelength Deviation 
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Peak Finding Recursion tt 
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MODEL OPM(x,H 0 ,p,Ac) :=Ho+ P (X -^c ) 2 Parabolic "near" the peak, at 

Lambda_c, Beta is negative 



Optical Power Meter Reading in dB 

Mj are the measured 
Optical Power values. 

J is the number of samples 
per spectrum, j = 0. . . J-l 

M spectra collected at 
T=4.5 second intervals 
and interpolated to 4 
seconds. 



H o is insertion loss 

P is inverse bandwidth squared parameter Samples of M can be collected 

X is center Wavelength individually, e.g. via laser 

heterodyne, and in any order. 


wfsA 
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Peak Finding Recursion 


Si 
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P(H,p,/r):=^ M. - H+p (x. -/c) 


H := 

/vw 


Performance Function: Find H ; Beta ; and Lambda_C 
(not initial values); After at least one iteration r 
REPLACE initial values H0 ; Beta_Q ; and Lambda_C0 
values with H ; Beta, Lambda C below; i.e. the 


P:= 





program is recursive and. optionally, iterative. 

z 

M J - p 0 )/.-/c 0 ) | 




j 


dBm 




J + 1 



BW = 3dB 


J 


BW" 

bandwidth 

z 

(M.-Ho^/.-^o) | 

dBm/pm A 2 

CT - 1 

Gaussian 
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^ 4 ln(10) p 

parameter 


Z 


^ „ \4 

/"j _AC °) 


- Ho - Po -fco) 2 ]-Po-(*. - /.c 0 )j 


/.c := AC t 


z 


2-p 0 2 (X. - /.C 0 ) 2 


M. -Ho-Po-p.j - AC 


of 


pm 


H -^Insertion Loss, p Bandwidth - 


A c -> Center Wavelength 


(the zero subscript denotes the value of the previous recursion) 
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A Barker 7 TheSIS Cycle for Fiber Bragg 

Gratings (FBG) 
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Temperature (blue) in deg C. AR3MA2 Model (red) 


T\ , ^ 0 ] ADlA/f AO ^R3 

MA2 

typical AKJMAz g2 g3 

AO A1 

Parameter Vector 1.5254s -0.2&380 -o.26sso 

0.81490 -0.81700 


T ch ip= 6.5 minutes 
T sampk= 4 seconds 


FBG Poles and Zeros 

introducing the "thermality zero" Z T 
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z-plane 



£ \Z- 


Context: T = 4 seconds/sample, Modified Barker 7 
“dither” sequence, Tspan = 6°C pp. 



Dominant FBG Poles and 
Thermality zero, Z T (Detail) 
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OS 0.9 1 11 12 
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Re(z) 


o 


I 2 3 4 5 6 7 8 9 10 11 12 13 14 


— RealPoleO 
-•-RealPolel 
-»-RealPole2 
-•-Zero Z T 



- 0.6 1 


Cycle # 


Stability of the FBG Thermality Zero, Z T 
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44 parts-per-million rms 
0.016% peak-to peak 




FBG Temperature Coefficient, 
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Cycle # 

NMA 

I Vi 

j zl 

NAR 

‘-I B . 


a 


<3% variation rms (600 kHz/°C) 

<9% variation peak-to peak (1.9 MHz/°C) 




Frequency Response Results 


j 


. s, 

Goddard Space 
Flight Center 


Frequency Response of FBG (mHz) 



Notes: 1) Mapped z plane to Laplace (s) domain via the bilinear transform 

2) There is a 4 milliHz peak (4 minute period) to be aware of. Thesis FBGs.avi 


AVI 
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Multivariate Figure of Merit 
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Example Target specification: 

A perfectly athermal DUT with a 
parameter vector comprising a specified 
center wavelength and bandwidth. 


Example DUT response: a DUT with a 
parameter vector comprising center 
wavelength, temperature coefficient, and 
bandwidth. 


Figure of Merit (FOM), i.e. the expected 
performance over the observed range of 
environmental conditions, e.g. 
temperature deviations. 

Note: o x represents a vector of parameters 
of the measured probability density 
function of temperature. 

Note: the Spec and DUT parameter 
vectors can include insertion loss (IL); i.e. 
one can penalize deviations from a (IL) 
specified value. 


i 
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( 2 Yl 

1 

-X 

■ exp 

2 . 

U'Ueq JJ 


wavelength here is relative to specification 


iJj(XAc-,at a cr ? I) := 


1 

-[X- (Xc + O-T)] 2 

eX P 

^i-TUCT 

2-cr 2 


measured wavelength deviation from spec 
and temperature coefficient 


1 span 


FOM*! Xc ? oc r cr r ctj i := 


Xspan 


(X,cr ) dX 


■ Xspan 


■PDFjT.Oj) dl 


- Tspan 


Given a specification, optimize FOM as a 
function of any DUT, thermal parameter, or 
combinations thereof. 



Optimum Temperature Coefficient of Wavelength 

Example: Asymmetric Temperature Regulation 



FOM Surface 

(independent 
variables = 
wavelength offset 
and bandwidth) 




Goddard Space 
Flight Center 


'A 



Summary 
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Points & Highlights: 

Demonstrated the use of spread spectrum temperature cycling to shorten test time. 

Tested models using the AICc and extracted parameter vectors suitable for control loop design. 
Introduced the idea of the “thermality zero”. 

Presented a recursive algorithm that extracts multiple parameters of passive DUTs. 

Showed 30 femtometer rms (la) precision in measuring the center wavelength of passive 
optical components (FBGs). 

Showed Model stability across multiple Barker 7 cycles: 2.8% rms cycle-to-cycle variation of 
thermal coefficient of FBGs (MHz/°C ) 


